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1. Introduction 

Formal random matrix models have been used for their interpretation as combinato- 
rial generating functions for discretized surfaces [4,5,24]. The hermitean one- matrix 
model counts surfaces made of polygons of only one color, whereas the hermitean 
two-matrix model counts surfaces made of polygons of two colors. In that respect, 
the 2-matrix model is more appropriate for the purpose of studying surfaces with 
non-uniform boundary conditions. At the continuum limit, the 2-matrix model gives 
access to "boundary operators" in conformal field theory [21]. 

Generating functions for surfaces with boundaries are obtained as random matrix 
expectation values. The expectation value of a product of I traces is the generating 
function for surfaces with I boundaries, the total power of matrices in each trace 
being the length of the corresponding boundary. If each trace contains only one type 
of matrix (different traces may contain different types of matrices), the expectation 
value is the generating function counting surfaces with uniform boundary conditions. 
Those non-mixed expectation values have been computed for finite n since the work 
of [16,25] and refined by [1]. 

Mixed correlation functions have been considered as a difficult problem for a 
long time and progress have been obtained only recently [2, 15]. Indeed, non-mixed 
expectation values can easily be written in terms of eigenvalues only (since the trace 
of a matrix is clearly related to its eigenvalues) , whereas mixed correlation functions 
cannot (TrM^M2' cannot be written in terms of eigenvalues of Mi and M2). 

The large N limit of the generating function of the bicolored disc (i.e. one 
boundary, two colors, i.e. < TrM^M2 >) has been known since [11,12,19]. The 
large limit of the generating function of the 4-colored disc (i.e. one boundary, 
4 colors, i.e. < TrM^M!^' M^" M^'" >) has been known since [8]. The all order 
expansion of correlation functions for the 1-matrix model has been obtained by a 
Feynman-graph representation in [9] and the generalization to non-mixed correlation 
functions of the 2-matrix model has been obtained in [10]. 

Recently, the method of integration over the unitary group of [15] has allowed to 
compute, for finite N, all mixed correlation functions of the 2-matrix model in terms 
of orthogonal polynomials. 

The question of computing mixed correlation functions in the large N limit is 
addressed in the present article. 

The answer is (not so) surprisingly related to classical results in integrable statis- 
tical models, i.e. the Bethe Ansatz. It has been known for a long time that random 
matrix models are integrable in some sense (Toda, KP, KdV, isomonodromic sys- 
tems,...), but the relationship with Yang-Baxter equations and Bethe Ansatz was 
rather indirect. The result presented in this article should give some new insight in 
that direction. We find that the /c-point functions can be expressed as the product 
of 2-point function, which is the underlying idea of the Bethe Ansatz. 
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Outline of the article: 

- section 1 is an introduction, 

- in section 2, we set definitions of the model and correlation functions, and we 
write the relevant loop equations, 

- in section 3, we introduce a Bethe Ansatz-like formula, and prove it in section 4, 

- in section 5, we solve the problem under a matrix form, 

- section 6 is dedicated to the special Gaussian case. 



2. The 2-matrix model, definitions and loop equations 
2.1 Partition function 

We are interested in the formal matrix integral: 

Z :^ [ ^ dM, dM2 e-^^ [^i(A^i)+^2{M2)+MiM2] (2-1) 

where Mi and M2 are N xN hermitean matrices and dMi (resp. dM2) is the product 
of Lebesgue measures of all independent real components of Mi (resp. M2). Vi{x) 
and V2{y) are complex polynomials of degree di + 1 and ^2 + 1) called "potentials". 
The formal matrix integral is defined as a formal power series in the coefficients of 
the potentials (see [5]), computed by the usual Feynman method: consider a local 
extremum of g-^'^ [^i(^i)+^2{M2)+MiM2]^ expand the non quadratic part as a 
power series and, for each term of the series, perform the Gaussian integration with 
the quadratic part. This method does not care about the convergence of the integral, 
or of the series, it makes sense only order by order and it is in that sense that it can 
be interpreted as the generating function of discrete surfaces. All quantities in that 
model have a well defined expansion [27]. 

The extrema of Vi{x) + V2{y) + xy are such that: 

Vi{x)^-y , Vl{y)^-x (2-2) 

there are did2 solutions (indeed V^'(— ^/(a;)) = —a;), which we note {xi,yj), 

7 = 1,..., did2. The extrema of T]:[Vi{Mi)+V2{M2)+MiM2\ can be chosen diagonal 

(up to a unitary transformation), with x/'s and y^'s on the diagonal: 

m times n2 times n^^ii^ times 



Ml = diag(xi, ...,xi,X2,...,X2,..., xa^di, ^didi) 

m times 712 times ^^^1^2 ^^^^^ 

M2 = diag(j/i, ...,yi,y2,...,y2,..., Vd.d^^ Vd^d^) (2-3) 

The extremum around which we perform the expansion is thus characterized by a 
set of filling fractions: 

did2 

1=1 
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To summarize, let us say that the formal matrix integral is defined for given 
potentials and filling fractions. 

The "one-cut" case is the one where one of the filling fractions is 1, and all the 
others vanish. This is the case where the Feynman expansion is performed in the 
vicinity of only one extremum. 

2.2 Enumeration of discrete surfaces 

It is well known that formal matrix integrals are generating functions for the enu- 
meration of discrete surfaces [4,5, 13, 18]. 

For instance, in the one-cut case (expansion near an extremum x, y), one has: 



lnZ = 

1 



^G#Aut(G)^^ \5) \S) \S) i^i=3 i/i lii=3 i/i 



where the summation is over all finite connected closed discrete surfaces made of 
polygons of two signs (+ and -). For such a surface (or graph) G, xiG) is its 
Euler characteristic, riiiG) is the number of i-gons carrying a + sign, fiiiG) is the 
number of i-gons carrying a — sign, n_|_+(G) is the number of edges separating two 

-I- polygons, n (G) is the number of edges separating two — polygons and n+_(G) 

is the number of edges separating two polygons of different signs. #Aut(G) is the 
number of automorphisms of G. 

The QiS, QiS and S are defined as follows: 



9k ■-- 



d''Vi{x) 



Qk 



8 :- g2~g2 - 1 (2-6) 



x=y 



Example of a discrete surface: 




(2-7) 



In the multicut case, i.e. with arbitrary filling fractions, matrix integrals can 
still be interpreted in terms of "foams" of surfaces, and we refer the reader to the 
appendix of [3] or to [14] for more details. 
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2.3 Enumeration of discrete surfaces with boundaries 

Similarly, given a sequence of signs Si, S2, ■ ■ ■ , Sk, Si G 1, 2, it is well known that the 
following quantity: 

(iVdlMj^ (2-8) 

is the generating function of discrete surfaces with one boundary of length k, whose 
signs of polygons on the edges are given by the sequence (si, . . . , s^). 

Example of a discrete surface with boundary (-|--|--|--|--|--| + + 

+ + ): 



(Tr (MfM|MfM|)) = E, 




(2-9) 



More generally, an expectation value of a product of n traces is the generating 
function for discrete surfaces with n boundaries. 

In this article, we are interested only in one boundary and to leading order in 
A^, i.e. surfaces with the topology of a disc. 

2.4 Master loop equation and algebraic curve 

Let us define: 

V^(x) := i- / T> \ , W^(7/):=^/Tr ^) (2-10) 

where the expectation values are formally computed as explained in the previous 
section, with the weight e-^^[^i(^i)+^2(M2)+MiM2]_ ^y^^.) ^^esp. W{y)) is defined 
as a formal power series in its large x (resp. large y) expansion, as well as in the 
expansion in the coefficients of the potentials. W{x) (resp. W{y)) is a generating 
function for surfaces with one uniform boundary, i.e. with only sign -|- (resp. sign 
— ) polygons touching the boundary by an edge: 



W{x) = 




(2-11) 
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We also define the following formal series: 

Y{x) := W{x) - V;{x) , X{y) := W{y) - V^{y) (2-12) 
In addition, we define: 



X - Ml 1/ - M2 x' -Ml/ 

-N'W{x'){U{x,y)-x-V^{y)) 

(2 - 15) 

y) := (l^/(x) + |/)(l^2'(l/) + + P{x, y) - 1 (2-16) 

Notice that U{x, y) and U{x, y; x') are polynomials of y (with degree at most ^2 — 1), 
P{x,y) is a polynomial of both variables of degree {di — l,d2 — 1) and E{x,y) is a 
polynomial of both x and y of degree (rfi + 1, ^2 + 1). 

It has been obtained in many articles [7, 11, 12,26], that: 

E{x,Y{x)) = ^U{x,Y{x),x) (2-17) 

To large N leading order that equation reduces to an algebraic equation for Y{x), 
called the "Master loop equation" [26]: 

E{x,Y{x)) = Q (2-18) 

(similarly, one also has E{X{y),y) — 0, which imphes y o X = Id, known as 
Matytsin's equation [23]). The coefficients of E{x,y), i.e. of P{x,y), are entirely 
determined by the conditions £4. ydx — 2mei for a choice of irreducible cycles on the 
algebraic curve. 

The properties of that algebraic equation have been studied in many works [7,20]. 
Here we assume that it is known. 

2.5 Correlation functions, definitions 

We define: 

W,{xi^yi^x,^...^x,^y,) := l^Tr A^;^^-^) (2-19) 
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Pol V^ivk) Wk{xi, yi, X2, . . . , Xk, Vk) 



Vk 

1 /^^ 1 1 1 V^iyk)-vm2) \ 

N\ ^ xi- Nhyi- M2'" Xk- Ml yk-M2 / 

(2 - 20) 

Pk{xi,yi,X2, . . . ,Xk,yk) 
:= Pol Pol Vl{xi) ^'(i/fc) Wkixi, yi, X2, . . . , Xk, Vk) 

^ / Vlixi) - VliMi) 1 1 V^iyk)-VKM2) \ 

N\ xi-Mi 1/1-M2 ■■■xfc-Mi yk-M2 / 

(2-21) 

Akixi, yi, X2, . . . , Xfc) := ^ ( Tr . . . L-^'(M2)\ (2-22) 

\ xi- Ml yi - M2 Xk - Ml I 

where Pol^^^ f{x) denotes the polynomial part at infinity of f{x) (i.e. the positive 
part in the Laurent series for x near infinity). 

The functions Wk are generating functions for discrete discs with all possible 
boundary conditions. One can recover any generating function of type eq. (|2-8|) by 
expanding into powers of the Xj's and ?/j's. 

For convenience, we prefer to consider the following functions: 

Wk{xi, yi, X2,..., Xk, yk) ■= Wk{xi, yi,X2,..., Xk, yk) + 4,i (2-23) 

Uk{xi,yi,X2,...,Xk,yk) ■= Ukixi,yi,X2, . . . ,Xk,yk) + + ^k) (2-24) 

and for k > 1: 

Pkixi,yi, X2,..., Xk, yk) ■= Pk{xi, yi, xa, . . . , Xk, yk) + Wk-i{x2, ...,Xk, yi) (2-25) 

For the smallest values of k, those expectation values can be found in the liter- 
ature to large N leading order: 

• it was found in [7, 11, 12]: 

Wi{x,y) = - f[^f , Uiix,y) = j^^^^ (2-26) 

[x - X{y)){y -Y{x)) {y-Y{x)) 

• it was found in the appendix C of [8] (there is a change of sign, because the 
action in [8] was ^-Ntr mMi)+V2{M2)-M^M2)y_ 

W2[xi, yi, X2, 2/2) = -, TT ^ (2-27) 

[xi - X2)[yi - 1/2) 

• For finite A^, it was found in [2], and with notations explained in [2]: 

Wi{x, y) = det M + ^N-i ^\^y^pJ ^N-i\ (2-28) 
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• For finite N, it was found in [15] how to compute any mixed correlation function 
in terms of determinants involving biorthogonal polynomials, with a formula very 



similar to eq. (|2-28 



Here, we shall find a formula for all WkS in the large N limit. 
2.6 Loop equations 

Loop equations are nothing but Schwinger-Dyson equations. They are obtained by 
writing that an integral is invariant under a change of variable, or alternatively by 
writing that the integral of a total derivative vanishes. 

The loop equation method is well known and explained in many works [7,26]. 
Here, we write for each change of variable the corresponding loop equation (we use 
a presentation similar to that of [7]). 

In all what follows we consider k > 1. 

• the change of variable: 6M2 = — Kj Kr ■ ■ ■ — imphes: 

k-1 

Ak{xi,...,Xk) = ^ W'j^xi, . . . . . . 

i=i 

^ xiWk-ijxi, yi,..., Vk-i) - XkWk-ijxk, 1/1, ... , ^fc-i) 

Xl Xk 

k-1 

= J2 ^ji^i^ • • • ' Vj) Wk-j{xk, Vj, . . . , Vk-i) 
i=i 

-Wk-i{xk,yi, . . .,yk-i) 

Wk-i{xi, 1/1, ... , yk-i) - Wk-i{xk, 1/1, ... , yk-i) 

+Xk 



(2 - 29) 



Xl - Xk 



the change of variable: SMi = — V? hn- ■ ■ ■ — j— ^afafc^ ^i'^^^^^ implies: 



2^1=2 11^7^ ^ k-]+l\Xj-, (/j, . . . , Xk, yk) 



- yk) Ukixi, ...,yk) 

i=2 

+v^{yk) 



^i=2 xi-Xj 

, Wh-i(xi,yi,...,yi:_i)-Wk-i(xk,yi,X2,...,yk-i) 

Xl-Xk 

+Ak{xi, ...,Xk)- Pk{xi, yi, X2,..., Xk, yk) (2-30) 

_ y-fc Wj^i(xi,yi,...,yj.i)-Wj-i{xj,yi,...,y^^i) jj , x 

— l^j=2 xi-Xj ^k-j+l\Xj, (/j, . . . , Xk, yk) 

+ EJ=J Wj{xi, yj) Wk-jixk, yj, yk-i) 
-Pk{xi,yi,X2, . . . ,Xk,yk) 



where we have used the loop equation eq. ( p^-29| ) for Ak{xi, . . . , Xk). 
• the change of variable: 5M2 = — Ki \t- ■ ■ ■ — in 

{X{yk) - Xl) Wk{xi, yi, X2, . . . , Xk, yk) 
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Wk-j{xj+i, ...,yk)- Wk^j{xj+i, ...,Xk, Vj) . 
— j{Xi, . . . , yj) 

-f/fc(xi, ...,yk) 

(2-31) 

2.7 Recursive determination of the correlation functions 



Theorem 2.1 The system of equations eq. { 2-3i ) and eq. ( 2-3]\ ) for all k has a 
unique solution. 

In other words, if we can find some functions Wk, Uk, Pk which obey eq. ( P-3(]| ) 
and eq. ( P-31| ) for aU /c, then they are the correlation functions we are seeking. 

proof: 

Wi, Ui and Pi have already been computed in the literature. 
Assume that we have computed Wj, Uj, Pj for all j < k. Let us show that 
eq. (|2-3CI|) and eq. ( |2-31| ) determine uniquely Wk, Uk and Pk- 



Let X^"\yk), a = 0, . . . ,di be the di + 1 solutions for x of E{x, yk) = 0. For 
every a = 0, . . . , (ii one has: 

Y{X'-^\yk))=yk (2-32) 
At xi = X(")(?/fc), eq. ( F30| ) reads: 

Pfc(X(")(?/fc),?/i,X2, . . ■,Xk,yk) 

= EU Uk-,^ii^„ ...,xk, yk) (2-33) 

+ Ej=! W,iX("\yk), . . . , %•) wi.ixk, . . . , yk-i) 

where all the quantities in the RHS are known from the recursion hypothesis. We 
thus know the value of Pk for di + 1 values of Xi. Since Pk is a polynomial in xi of 
degree at most di — 1, we can determine Pk by the interpolation formula: 

. ^^ Pk{xu...,yk) 

^^^-""^y^y^ Eix.,yk) 
^ (X(")(|/fc) - Xjyk)) PkiX(-\yk), ...,yk) 
(xi-XH(|/fc))E.(XH(yfc),y,) 

(2 - 34) 



where Xk = X{yk) denotes X'^^^yk). Once Pk is known, eq. (|2-30|) allows to compute 
Uk-i and eq. (|2-31|) allows to compute Wk- • 

3. A Bethe Ansatz-like formula for correlation functions 

Thus, the loop equations determine Wk uniquely, i.e., if we can find Wk, Uk = 
PoW Vi{yk)Wk and Pk = Wk-i + Pol,, Vi{x{)Uk which satisfy eq. (ElD.eq. ( F3q) , 
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it means that we have the right solution. We can thus make an Ansatz for Wk, and 
check that it satisfies the loop equations above. 

Our Ansatz is similar to the Bethe Ansatz [17]: 

k 

Wkixi,yi,...,Xk,yk)= Cl^\xi,yi,---,Xk,yk) Y[Wi{xi,y^(^i)) (3-1) 

where the coefficients Ca are rational fractions of the Xj's and i/j's, with at most 
simple poles at coinciding points, and independent of the potentials. 

We call eq. ( ^-1| ) a Bethe Ansatz-like formula, because it is very similar to the 
solution initially found by Bethe for the 1-dimensional spin chain, and then for the 
(5-interacting bosons. 

If we assume that eq. ( p-l| ) satisfies eq. ( p-31| ), we can in particular take the 
residue of eq. ( P-31| ) at yk Y{xi) for some /. That implies the following relationship 
among the coefficients C^'^^^s: 

(xa-^ik) - xi) Cl'^^xi, yi,..., Xk, yk) 

^ ^ , Ci^Kx,,...,y,)C(;-^\x,^u...,yk) 

j=irGS(i,...,j)pGS{i+i,...,fc) yj 

(3-2) 

Beside, since Wk is the expectation value of a trace, the C^'^^'s must be cyclically 
invariant: 

na;i,2/i,a;2,l/2, • • -^Xk.yk) = Cl''\x2,y2, ■ ■ . ,Xk,yk,xi,yi) (3-3) 
and, since Wk should have no poles at coinciding points yk = yj one should have: 

Res Ha^i, 1/1, X2, 1/2, • • • , Xk, y'k)dy'k = Res CfJ;L^(a;i, yi, X2, y2, ■ ■ ■ , Xk, y'k)dy'k 

(3-4) 

With C^^^ = 1, it is clear that the set of equations eq. (P-2|), eq. (|3-3|) , eq. ( P-4|) , have 
at most a unique solution. We prove in the next section that the solution exists, and 
thus, eq. ( |3-2| ), eq. ( |3-3| ), eq. ( |3-4| ) determine C^^^ uniquely. The C^'^^'s are explicitly 
computed in section ^. 

Then in section ITT, we prove that: 

Theorem 3.1 // the C^^^^'s are rational functions defined by eq. ( \3-3i ), eq. I\3-^ ), 
eq. ( \3-4\) , then the functions Wk 's defined by the RHS of eq. ( \3-]\ ), the functions 
Ukixi, . ..,yk) ■■= Fo\{V2{yk))Wkixi, ...,yk), and the functions Pk{xi, ■ ■■,yk) ■ = 
Po\(yl{xi))JJk{xi, ...,yk) + Wk-i{x2, ...,Xk, yi), satisfy eq. §^ and eq. (f^. 



As a corollary, using theorem p.l|, we have: 
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Theorem 3.2 




(3-5) 



The derivation of theorem |3.1| is quite technical and is presented in section 



4. Amplitudes of permutations 

In this section, we compute the amplitudes C^'^^ explicitly. 

Eq. ( p-2|) , eq. (p-3|) , eq. (1^1) ^^(^ initial condition C^^\xi,yi) = 1 clearly define 
at most a unique function C^^\xi, . . . ,yk). In this section, we build the solution 
explicitly, and then, we prove that the function we have constructed indeed satisfies 
eq. (H), eq. O, eq. (0). 

Below we prove that eq. ( |3-2D implies that C^'^-* vanishes for non planar permuta- 
tions (see Definition 4.1 below) and, for planar permutations, C^^'^ is the product of 

(k) 

Cj\j corresponding to faces. We are thus led to introduce the following definitions: 

4.1 Some definitions: planar permutations 

Let S be the shift permutation: 

^:=shift = (1,2,...,A;-1,A;) , i.e. S{€)=i + 1 (4-1) 
Definition 4.1 A permutation a ET^k is called planar if 

?^cycles(cr) + 'T'cycles (^'o") =k + l (4-2) 

where ncycies(c") the number of irreducible cycles composing the permutation a. 
Let Sfc C Sfc be the set of planar permutations of rank k. 

Eq. is equivalent to saying that if one draws the points Xi,yi, . . . ,Xk,yk 



on a circle, and draws a line between each pair {xj,ycr{j)), the lines don't intersect. 
The cycles of a and the cycles of Sa correspond to the faces (i.e. the connected 
components) of that partition of the disc. 

Each planar permutation can also be represented as an arch system, and thus, the 
number of possible planar permutations is given by the Catalan number Cat {k): 

— 

Card(E,) = Cat {k) = ^, (4-3) 
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X ^0(1) X„„,,, 

a(l) _ Sci(l) 




o o-cycle 
So-cycle 



(4-4) 



Example of a planar permutation and its faces. 
4.2 Face amplitudes 

Definition 4.2 For any k > 1, we define a rational function of xi, . . . ,yk: 

F^''\xi,yi,X2,...,Xk,yk) (4-5) 



by the recursion formula: 

F«(xi,i/i) : 
F^''\xi,yi,...,Xk,yk) ■ 



F^^\xi,yi, . . .,Xj,yj)F^'' ^\xj+i,yj+i, . . .,Xk,yk) 



(4-6) 



Lemma 4.1 For any k > 1 



E 



{xk - xi){yk - yj) 



FW(xi,yi,...,Xfc,i/,) = 0(i/^'^) 



(4-7) 



when yk ^ oo. 



proof: 

Let us prove it by induction on k. It is true for k=l. Let k be larger or equal to 2 
and assume that this is true for all F^^'^ with j < k. Then eq. ( |4-6| ) straightforwardly 
gives the same behaviour for F^^\ • 



Lemma 4.2 F^'^'^ has cyclic invariance, i.e. 

F^''\x2,y2, ■ ■ .,Xk,yk,xi,yi) = F^^\xi,yi, . . .,Xk,yk) 



(4-8) 



proof: 
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We prove it by recursion. It is clearly true for k = 1 and k = 2 since 

X2, 1/2) = ■ For k > 3, assume that it is true for all F^^'^ with 

j < k. One has: 



{xi -X2){yi -yj) 



i=2 

j^ F(^-^\x2,...,x„y,) 
P2 {xi - X2){yi ~ yj) 

yi, Xj+i, ^j+i, ...,Xi, ■ ■ ■ , Xfc, ^fc) 

(a;fc-a;i)(yfc-yi) 
^ ^ F(-^'-^)(x2, ■ ■ ■ , Xj, yj) F(^) (xi, yi)F^''-^\xj+i, yj+i, ...,Xk, yk) 
Pi {xi~ X2){yi-yj) {xk - xi){yk - yi) 

gF('=-')(a;;+i,|/;+i,...,Xfc,|/fc) 



i=3 



(xfc - xi){yk - yi) 

^ F(j-^)(x2,...,X j, . . . , X,, y;, Xi, 

i=2 



(Xi -X2)(l/l 

^1 F(^-')(x;+i, ■ ■ ■ , Xfc, yk)F^^\x2, • ■ • , X;, |/;, Xi, 

F('^(xi, 7/1, X2, ■ ■ ■ , xi, yi)F^^-^\xi+i, yi+i, ...,Xk, Vk) 
{xk - xi){yk - yi) 



F''^\xi,yi, . . .,Xk,yk) 



(4-9) 



Lemma 4.3 For k >2, F^^'^ has simple poles in yk'- 

k-l 

F^''\xi,yi,...,Xk,yk) = Y1 Res F^''\xi,yi, ... ,Xk,yk)dy'k (4-10) 

yk - yi y'k^vi 

proof: 

It is clearly true for k = 2. We prove it by induction on k. Assume that it is 
true up to A; — 1. Using the recursion hypothesis, one can see that each term in the 
RHS of eq. (13) has at most a simple pole at yk = yi and one can write eq. ( [4-10| ) 
with the use of Lemma [4.1| . Thus the recursion hypothesis is true for k. • 

4.3 The amplitudes 

Definition 4.3 Then, for any k > 1, and for any permutation a G S^, we define 
C^'')(xi, yi, X2, . . . , Xk, yk) a rational function of Xi, . . . , yk, by: 
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• C^''\xi, ui, X2. . . . .Xk,yk) '■= if a is not planar, and 

• if a is planar, we decompose a and Sa into their product of cycles: 

a = aia2...ai , Sa = aia2 ■ ■ ■ ai (4-11) 

such that: 

^3 = (Hu h,2, • • • , , (^{'ij,m) = ij,m+i (4-12) 

= H2, , (^i'ij,m) = Hm+i - 1 (4-13) 



I 

Ci''\xi,yi, X2,..., Xk, yk) := H ^^''^(^i^.i' 2/ij.2, a;i,,2, |/i^,3, • • • , Xi^^, , |/i .J 

(4 - 14) 

i.e. C^^ is the product of F's of each connected component of the disc partitioned by 
a. 



4.4 Examples 

In particular with a — Id, we have: 

Cidi^i, yi, X2,..., Xk, yk) = F'^''\xi, yi, X2, y2, ■ ■ ■ , Xk, yk) (4-15) 
and with a — S^^: 

C^s-i(xi, yi, X2,..., Xk, yk) = F^^\xk, yk-i, Xk-i, ■■■,yi, xi, yk) (4-16) 
An example with k — 12 

Let us consider an example of cr e E12 defined as follow: a{l) — 3, a{2) = 1, 
ct(3) = 2, ct(4) = 7, ct(5) = 6, ct(6) = 5, a{7) = 4, a{8) = 8, a{9) = 11, (7(10) = 10, 
(7(11) = 12 and (7(12) = 9: 

(7 = (1,3,2)(4,7)(5,6)(8)(9, 12,11)(10) 
Sa = (1, 4, 8, 9)(2)(3)(5, 7)(6)(10, 11)(12) (4-17) 
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The corresponding arch system is: 




(4-18) 



where dark faces (resp. white faces) correspond to the cycles of a (resp 
that permutation C^^^^ is worth: 

^^^Ha^i' y3, 3^3, y2, X2, yi)F(2)(a;5, ye, a^e, ys)^^^^ (a^g, yi2, a;i2, ^11,2:11,^9) 
xF(2)(x4, y7, 2^7, y4)i^^^^ ys, 3^4, y?, 2^8, ys, a^g, yi2)i^^^^ (a^s, ye, a;?, y4) 
xF(2)(xio,yio,a;ii,y9) 



. So). For 
(4-19) 



Example A; < 3 

C(^) = 1 
1 



C, 



(2^2 - a;i)(2/2 - 2/1) 

(2) ^ \ 

(''^ (^2-Xi)(|/i-2/2) 



(xi - X3) Vys - yi (x2 - 3:^3) (ya - Vn) ya - 2/2 (a^^i - a:;2)(y2 - Vi) , 



(3) 

(1)(23) 



11 1 

(xi - X2) 2/3 - 2/1 (a^2 - 2:3) (2/2 - 2/3) 



^(3) ^ 1 1 

{12){3) _ _ (3,^ _ 3^^) ^y^ _ ^ 

C<(3) ^ ^ ^ + ^ ^ 

( xi - X2) (x2 - X3) (yi - ^2) ya - yi (a^i - a^a) (a^i - 2:2) (y2 - yi) ya - y2 

^(3) ^ I 

^^^^^^^ ya - yi (xi - X3){x2 - X3){yi - 2/2) 

C(?32) = (4-28) 
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Example: Rainbows 

The rainbow is the permutation 

a{j) = k + l-j (4-29) 

if k is even: 

ci^) = L 

Uiilixk+i-i - Xi){yi - Vk+i-i) I\i=i~^{xk+i-i - Xi+i){yi - yk^i) 

(4 - 30) 
if k is odd: 

I\fji^''^{xk+i-i - Xi){xk+i-i - Xi+i){yi - yk+i-i){yi - yk-i) 

(4-31) 
4.5 Properties of 

Lemma 4.4 The C^^'' 's are cyclically invariant: 

Cl''\xi,yi,X2,y2, ■ ■ ■,Xk,yk) = Cl''\x2,y2, ■ ■ . ,Xk,yk,xi,yi) (4-32) 

proof: 

It follows from Lemma J4.2|. • 



Lemma 4.5 The C^'^^ 's have at most simple poles in all their variables and are such 
that: 



Ci^\xi,yi,X2,y2,---,Xk,yk) = J2 P-^^ Cl''\xi,yi,X2,y2, ■ ■ ■ ,Xk,y'k)dyk 

yk - yj y'k^yj 

(4-33) 

proof: 

If a is planar, the pair {ykiUj) can appear in at most one factor of eq. ( |4-14| ), 
and the results follows from Lemma 14.31. • 



Theorem 4.1 The C^^'^ 's, with a(l) ^ k, satisfy the recursion formula eq. (\3-2j : 



I t:^ ^ _ . 1 



C('^) = E E 5^,r,^^Ci^^Cf-'^ (4-34) 

Xa ^{k) ^1 j = l reS(l,...,j) pgS(i+l,..,fc) yj 

proof: 

Let TT be the cycle of length / -|- 1 > 2 of Sa which contains Xi and yk'- 

TT = (xi ^ yi^-i Xi^ yi2-i ^ . . . ^ ^ yk xi) (4-35) 
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ij+i ■.= a{ij) + l , io:=l (4-36) 

Planarity implies that: 

io < ii < i2 < ■ ■ ■ <ii (4-37) 
Since a is planar, there exists a unique way of factorizing a as: 

I 

I [ , (Tj G E(zj, . . . , ij+i - 1) (4-38) 

j=0 

From the definition of C^*^-* we have: 

^(fc) = F('+i)(xi,i/i,_i,Xi,,y,,_i, . . . ,i/i,_i,x,,,?/fc) nC^'+'^^'Ha;*,, • • • ,l/i,+i-i) 

(4-39) 

and using eq. ( [4-6| ), we have: 

= f('+i)(xi, Xii, . . . , a;^;, i/fc) Uj ciT'~''\xi^, yi^+,-i) 



notice that ii = a^^{k), and note: 



(4-40) 



m—l I 

Tm ■= n ' •= n (4-41) 

j=l j=m 

We have: 

G S(l, . . . , - 1) , p„ G S(i„, ...,k) (4-42) 



Eq. (1^) gives: 



1 i 1 

^(fc) _ y QMMk-ir.) /4_43N 

It is clear, from the planarity condition that if there exists some j and r and p, such 
that: 

a = Tp , rGE(l,...,j) , p G S(j + 1, . . . , fc) (4-44) 

then, one must have j = im, t = Tm and p = Pm for some m. • 

Lemma 4.6 For any transposition {k,j) (with k ^ j), we have: 

Res Cl''\xi, 2/1, X2, 2/2, • • • , 2/^)^?/^ = - Res Cf^L^(a;i, yi, Xs, 2/2, • • • , a;^, 2/^)^2/^ 

(4-45) 
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proof: 

It is trivial if Uj and yk are not in the same face, because both sides vanish: the 
LHS has no pole and the RHS is a non-planar permutation. The case where Uj and 
Uk belong to the same face reduces to proving the Lemma for a = Id. 

For a = Id, we prove it by recursion on A;. It clearly works for k = 1 and k = 2. 
Assume that it works up to A; — 1. 

From the definition eq. (|4.3| ) we have: 



yi, X2, 1/2, • • • , Xk, Vk) 
= F^^\xi,yi, . . . ,Xj,yk)F^^~^\xj+i,yj+i, . . . ,yk-i, Xk,yj)F^'^\xj,yk, Xk^Vj) 
F^^\xi,yi, ■ ■ ■ , Xj, yk)F^^-^\xj+i, yj+i, . . . , yk-i, Xk, yj) 
{xj - Xk)iyk - yj) 

(4 - 46) 



and thus: 
P-es dyjxi,yi,. . . ,Xk,yk) 



[k) , ^ _ ■ ■ .,Xj,yj)F^^ . . .,Xk,yj] 



(4-47) 

On the LHS, we have from eq. ( [4-6|) : 



Res F^''\xi,yi,X2,y2,---,Xk,y'k) 

yi, X2, . . . , xi, yi)F^^-^\xi+i, yi+i, ...,Xk, y'k) 

= ties } -. -— r 

y'k^yjj^i [xk- xi){y^-yi) 

^ F^^\xi,yi,X2, ■ ■ ■ , Xj, yj)F(^~J')(xj+i, yj+i, ■ ■ ■ , x^, yj) 

{Xk - Xi) 

^■^ F^^\xi,yi,X2,...,xuyi) ^f^n 
+ 2^ 7 77 ^ Res F^ '{xi+i, yi+i, ...,Xk, yk) 

[Xk - Xi){yj - yi) y'^^yj 

(4 - 48) 

Then, from the recursion hypothesis, and from eq. ( [4-4 7| ) we have: 

Res F^^\xi,yi,X2,y2, ■ ■ ■ ,Xk,y'k) 

F^^\xi,yi, X2,..., Xj,yj)F^^-^\xj+i, yj+i, ...,Xk, yj) 

{Xk - Xi) 

F^^Kxi,yi,X2,...,xi,yi) Mk-i)^ /^ 
F(^)(xi, 1/1, X2, ■ ■ ■ , Xj,yj)F'^''~^\xj+i, yj+i, ...,Xk, yj) 

{Xk - Xi) 

F^^\xi,yi,X2,...,xi,yi) 
ir[ (xk- xi){yj -yi) 

F^^-^\xi+uyi+u ■ ■ ■ ,Xj,yj)F^''-^\xj+i,yj+i, . . . ,yk-i, Xk,yj) 

i^Xj Xk) 
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(4 - 49) 

In the last line, we use again eq. (|4-6|) and get: 



Res F^''\xi,yi, X2, 2/2, • • • , Xk, y'k) 
^ F^^\xi,yi,X2, . ■ .,Xj,yj)F^''--^\xj+i,yj+i, . . .,Xk,yj) 

(Xk - Xi) 

_r _ ^. F'^^Kxi,yi,X2, ■ ■ ■,Xj,yj) F^''-^\xj+i,yj+i, . . ■ ,yk-i, Xk,yj] 

^ ^ {Xk- Xi) {Xj - Xk) 

(4 - 50) 



One can then see that the sum of eq. ( [4-47| ) and eq. ([4-50|) vanishes, and the recursion 
hypothesis is proven for k. • 

Lemma 4.7 

V/Ol, 5] Cf)(xi,|/i,X2,2/2,...,a;fc,i/fc) = (4-51) 

proof: 

This expression is a rational function of all its variables. Consider the poles at 
yk = yj, and write r = (/c, j). One can split the symmetric group into its two 
conjugacy classes wrt the subgroup generated by r: = [Id] © [r]. In other words: 

C„{xi,yi,X2,y2,---,Xk,yk) 

= C„{xi,yi,X2,y2, ■ ■ ■ ,Xk,yk) + Cra{xi,yi,X2,y2, ■ ■ ■ ,Xk,yk) 

(4 - 52) 



From Lemma |4.6| , the terms in the RHS have no pole at yk = yj- Similarly, using 
cyclicity and doing the same for the x's, we prove the lemma. • 



The Lemmas we have just proven, are sufficient to prove the main theorem 3.2 



This is done in section 4.7 



4.6 Computation of the rational functions F^''\xi,yi, . . . ,Xk,yk) 

Although the exact computation of the F^'^^'s is not necessary for proving theorem 



37^ , we do it for completeness. In this section we give an explicit (and non-recursive) 
formula for the F^'^^'s. 

A practical way of computing these formulas is described in Appendix A. 

Definition 4.4 To every permutation a G T^k-i, we associate a weight f^r computed 
as follows: 

I Ifi I Ifi /i 

= n n 9i„,l,ir^,J,^n,J + l II ll S'in j Jn.l -f^^n j ) ll 9iij,k,a(il_j) (^'SS) 

n=lj=2 n=2j=2 j=l 
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where gt^hj is defined as follows: 



1 



1 



gi,h,j ■= (4-54) 

Xh - Xi yh- Vj 

and a and Sa are decomposed into their product of cycles as in Definition K-Sj . 

Theorem 4.2 F^''\xi, yi, . . . , Xk, yk) is obtained as the sum of the weights fa 's over 
all a G 

F^^\xi,yi,...,Xk,yk) = J2 /- (4-55) 

proof: 

First of all, let us interpret diagrammatically the recursion relation eq. ( [4-6| ) 
defining the F's: 





(4-56) 

Actually, this recursion relation is nothing else but a rule for cutting a graph 
along the dashed line into two smaller ones. The weight of a graph is then obtained 
as the sum over all the possible ways of cutting it in two. 

Notice that F^''^ is the sum of Cat {k — 1) different terms. 

Let us now explicit this bijection with the graphs with k — 1 arches. In order 
to compute one of the terms composing F^''\ one has to cut it with the help of the 
recursion relation until one obtains only graphs with one arch. That is to say that 
one cuts it — 1 times along non intersecting lines (corresponding to the dashed one 
in the recursion relation). If one draws these cutting lines on the circle, one obtains 
a graph with k — 1 arches dual of the original one. Thus every way of cutting a graph 
with k — 1 arches is associated to a planar permutation a G S(i . . .k — 1). Let us 
now prove that the term obtained by this cutting is equal to fa- 

For the sake of simplicity, one denotes the identity graph of {xj, yj, . . . ,Xk, yk) 
by circle {j,j + 1, . . . , k). In these conditions, the recursion relation reads: 



;i,2,...,fc) = ^(7i,,,,(i,...,j)(j + i,...,fc) 



(4-57) 



Let 0" be a permutation of (1, . . . , A; — 1). Cut it along the line going from the 
boundary {xi,yk) to {ya(i),xsa{i))- It results from this operation the factor gi,k,a{i) 
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and the circles (1, . . . , and {Sa{l), . . . ,k): 

{l,...,k) . . . , ail))iSa{l), ...,k) (4-58) 

Then cutting the circle {Sa{l), . . . ,k) along {yk,xsa{i)) ^ {yaSa{i), xsaSa{i)) 
gives: 

{Sa{l), ...,k)^^ 9sa(i),k,aSa(i){Sa{l), aSa{l)){SaSa{l), ...,k) (4-59) 

One pursues this procedure step by step by always cutting the circle containing 
k. Using the former notations, this reads: 

{l,...,k)^^l[ gi,„k,a(l,^(ii,j. '^{iij)) (4-60) 

So one has computed the weight associated to the first Sa - cycle. The remain- 
ing circles correspond to cr-cycles. Let us compute their weight by considering for 
example (ii,i, . . . , (7(ii,i)) = (ii,i, . . . , ^1,2). 

The cut along the line {xi^^-,, yi^^^) ^ (z/n.s^ 

il,2) 5n,i,n,2,n,3(^l,l' • • • ' ^l,3)('5(ii,3), • • • , ^1,2) (4-61) 
Keeping on cutting the circle containing ii 1 at every step gives: 

h 

■ ■ • ,^1,2) n 9iiM,iij+i{S{tij+i), . . (4-62) 

One can notice that the remaining circles in the RHS correspond to cycles of 
Sa whose contribution has not been taken into account yet. One can then compute 
their values by following a procedure similar to the one used for the first Sa-cycle. 

One can then recursively cut the circles so that one finally obtains only circles 
containing only one element. This recursion is performed by alternatively processing 
on (T-cycles and S'a-cycles. 

Thus, one straightforwardly finds: 

I In I In ^1 

(1, . . . , A;) Yi n 9in,l,in,j,in,j+l 11 IT 9i„,j ,in,l,a{in,j) II 9i^,j ,k,a{hj) = U (4-63) 

n=lj=2 n=2j=2 j=l 

And then: 

F(^)= ^ U (4-64) 
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Example: Let us compute the weight associated to the permutation a G S12 
introduced earher. Starting from the circle (1, . . . , 13), one will proceed step by step 
the following cutting: 




(4-65) 



The first step consists in cutting along the 5"i cycle. The dashed lines show where 
one cuts the circles. Note that one do not represent the circles of unite length. The 
associated weight is 5(1,13,3 54,13,7 58,13,8 5'9,i3,i2- 

The second step consists in cutting along the remaining a cycles. One associates 
the weight 51,3,2 51,2,1 x 54,7,4 x 59,12,11 59,11,9 to this step. 

The weights associated to the two last cuttings are 5(5,7,6 x 5(10,11,10 and 5(5,6,5- 




(4-66) 



Finally, the weight of this planar permutation is then: 

fa — 9'1,13,3 9'4,13,7 5'8,13,8 5'9,13,12 5'l,3,2 5'l,2,l 5'4,7,4 5'9,12,11 5'9,11,9 5'5,7,6 5'lO,ll,10 9'5,6,5 

(4-67) 
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4.7 Proof of the main theorem 



We now prove that the function W defined by the RHS of eq. ( |3-1| ) and the functions 
U and P defined in theorem satisfy the system of equations eq. (|2-30|) and eq. 



Proof of theorem |3.1| : 

Using eq. (|4-33|) , one has: 



Ukixi, ...,yk) = Pol V^ivk) Wkixi,yi, ...,Xk, Vk) 

Vk 

k 

= Pol^'(?/fc) Cl''\xi,yi,...,Xk,yk) Y[Wi{xi,y„^i)) 

crGSfc i=l 

k-1 

= H Z] CaK^uVu- ■ ■^Xk,yk)dy'k Y[Wi{x^~m),yi) 

o-eEfe jVfc i=i 

y2{ykWi{x^-^k),yk) 
y" yt - yj 

k-l 

= J2 ^es Ci''\xi,yi,. . .,Xk,y'k)dy[ Y[Wi{x^-i^i),yi) 

o-eSfc jj^k ^k^yj i=i 

p^^ {W{yk) - X{yk))Wi{x^-^k).yk) 

y^ yk - yj 

= - II II Res (:7f)(xi,?/i,...,Xfc,?/fc)c??/fc X{Wi{x„-n^i),yi) 

o-eSfc jVfc yk-^yj i=i 

^^^ X{yk)Wi{x^-i(k),yk) 

y'' Vk - Vj 

k-l 

= H H Cl''\xi,yi,...,Xk,y'k)dy[ YlWi{x^-^i),yi) 

o-GEfe j^k Vk^Vj i=l 

p^^ {x^-H^k) - X{yk))Wi{x„-i(^k)^yk) - {.x„~i(^k) - X{yj))Wi{x„~i(^k)^yj) 

yk - yj 

k-l 

= H Res C^^^(xi,?/i,...,Xfc,?/^)rf?/fc X{Wi{x„-i(i),yi) 

crgEfc j^k yk^yj i=i 

{xg-^k) - X{yk))Wi{x„-i(^k),yk) - {xa-^k) - X{yj))Wi{x„-i{k),yj) 

Vk - Vj 

(4 - 68) 

Indeed, using eq. ( |2-26| ), one sees that the last expression is a polynomial in y^. 

ixg-i(^k) - X{yk))Wi{xg-^k),yk) - (a;<x-i(fc) - X{yj))Wi{xg-^k),yj) 

yk - yj 

-^(^^-i(fc)-yfc) _ ^(^a-i(k)'yj) 
yk - yj 

(4 - 69) 
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We have to check eq. ( |2-3 1|) , i.e. that A = with: 
^ j2 Wk-j{xj+i, ...,yk)- Wk-jjxj+i, ...,Xk, Vj) — 

-Uk{xi, . . ■,yk)^ 

+ - X{yk)) Wk{xi, yi, X2, . . . , Xk, yk) 

(4 - 70) 
We have: 

^ = E E E ^^^'^^(^^•+^' ■■■^y^^ 

j^k TT r yk — yj 

j k 

xY[Wi{Xr-^i),y,) Yl Wi{x^-^i),yi) 

i=l i=j+l 

-l^l^l^ 77—, [x^-Hk) , 

^\{Wi{Xr-l{i),yi) n W^l(a;7r-l(i)5Z/i) 

-2^2^ r^— W^{x^-i(^k), Vk) 

a j-tk yk Uj 

Res ) nVi(a:,-i(,),y,) 

fe-i 



X 



X 



-{xi-X{yk)) Y,C^^\xi, . . . ,yk)Wi{x„-^k),yk) H Wi{x„-i(i),yi) 



(4-71) 



i=l 



Using eq. ( |4-33|) in the last hne, adding it to the 4th hne, and using eq. ( [4-33| ) again, 
we get: 

^ = E E E ' " ' • • • ' 

3 k 

^X{Wi{xr-i{i),yi) n W^i(a;7r-i{i),Z/i) 

- E E E ^r\^i,---,y3)c^r^\^,^i.---.y3) w,{x^-.^k),y,) 

j^k -n- T Vk yj 
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X 

(4 - 72) 



i k-l 

i=l i=j+l 

+ 1^1^ Wi{x^-^k), Vj) 

k—\ 

^ i=i 



Using eq. ( |3-2| ) in the second line, exactly cancels the first line, and thus we get: 

A = - E E E • • • ■ -'"^ y,) 

j=ik T T Vk Uj 

j k-l 

i=l *=J+1 

+ 2^2^ Wi{x^-i^k), Vj) 

a yk yj 

k—l 

Res nW^i(^--Mi)>l/0 



X 

(4 - 73) 



i=l 



which is a rational fraction in y^ with poles at y^ = yj for some j. From Lemma 
A as defined in eq. ( |4-70| ) cannot have poles at yk = yj, thus A = 0. 

Now, we have to check eq. ( |2-30|) Using eq. ( |4-33|) , one has: 

Pfe(xi, ...,yk)- Wk-iix2, ■ . ■,Xk,yi) 
= Po\Vl{xi)Uk{xi,yi,...,Xk,yk) 

k—l 

= PolF(Xi) E Res llWl{x.-H^),y^) 

o-eSfc j^k yk-^Vi i=i 

Ui{xa-i{k),yk) - Uiix^-^k),yj) 
yk - yj 

k—l 

= Poir(xi) E Res c^'^) nw^i(^--W'2/0 

Uiixi,yk) - Ui{xi,yj) 

Vk ~ Vj 

+ PolF(xi) Y E Res l^i(a;i,y.(i)) H W^{x^-^^^,y■, 
^i(a^<7-i(fc),yfc) - Ui{Xa-i{k),yj) 

yk - yj 

= - E EE R^i Res nw^i(^-Ho>2/0 

c7eSfe,(T(l)=fc jT^fc ^7^1 ^ ' yk^V] i=l 
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E{xi,yk) - E{xi,yj) - E{xi, yk) + E{xi, yj) 
(xi — xi){yk — yj) 

+ E EE R^i Resc(^) n 

Xi — Xi 

Uiix^-Hk),yk) -Ui{x^-nk),yj) 

Vk - Vj 

(4 - 74) 

In order to satisfy eq. ( P-3CI|) , we must prove that 5 = 0, where: 

^ j2 Wi-i{xi,yi,...,yi_i) -Wi_i{xi,yi,X2,...,yi-i) 

1=2 ^1 ~ 

xUk-i+i{xi,yi, . . . ,Xk,yk) 

k-l ^ 

+ J2 Wii^u ■■■^yi) Wk-iixk, yi,..., yk-i) 
1=1 

-Pk{xi,yi,X2, ■ ■ ■,Xk,yk) 
-(F(xi) - yk) Ukixi, ...,yk) 

(4 - 75) 

One does it in a way very similar to the previous one, i.e. first prove, using 
eq. (13)' ^^^^ -B is a rational fraction of Xi, with poles at Xi = xi. But B can have 
no pole at Xi = xi, so,B = 0. 



5. Matrix form of correlation functions 

So far, we have computed mixed correlation functions with only one trace, i.e. the 
generating function of connected discrete surfaces with one boundary. In this section, 
we generalize this theory to the computation of generating functions of non-connected 
discrete surfaces with any number of boundaries. In order to derive those correlation 
functions, a matrix approach of the problem, similar to the one developed in [15], is 
used. 

Definition 5.1 Let k be a positive integer. Let ir and it' be two permutations ofY^k 
and decompose tt'^^-k into the product of its irreducible cycles: 

n'-\ = P,P2...Pn (5-1) 
Each cycle Pi o/tt'^^tt, of length pi, is denoted : 

J-) / • TT ' 7r'~^ ' TV ' TV^~^ Tv'~^ ■ TV ' 7v'~^ ' \ 

m V"m^l ^ Jm,l '"^ ^m,2 ^ Jm,2 • • • '"^ '^m^pm ^ Jm,pm. '"^ ^m,ly 

(5-2) 



26 



For any (xi, yi, X2, 2/2, • • • , Xk, yu) G C, we define 



-| Pm 



\m=l 



(5-3) 

which is a k\ X k\ matrix. 

Let us now generalize the notion of planarity of a permutation. 

Definition 5.2 Lei k be a positive integer. Let n and tt ' be two permutations ofJ^k. 
A permutation a & is said to be planar wrt (tt, tt') if 

ncycies(7r"V) + racycies(7r'"V) = k + ncydcs (tt'" V) (5-4) 

Let S^j^''^ ■* C Sfc 6e the set of permutations planar wrt (vr,7r'). 

Graphically, if one draws the sets of points (xj^ yj^^^, Xi^^^, 2, . . . , a:^^^^ , Vji^p^), 

n circles and link each pair {xj,ycr(j)) by a line, these lines do not intersect nor go 
from one circle to another. 

Remark 5.1 One can straightforwardly see two properties of these sets: 

• This relation of planarity wrt to (vTjVr') is symmetric in vr and tt' , that is to say: 

• The planarity defined in eq. ( [4.1| ) corresponds to tt = Id and tt' = S^^: 

Efc = Sf'^") (5-6) 

Directly from these definitions and the preceding results comes the following 
theorem computing any generating function of discrete surface with boundaries. 

Theorem 5.1 




(5-7) 

where is the k\ x k\-matrix defined by: 

• C^^^,{xi,yi,X2,y2, ■ ■ ■,Xk,yk) := if a is not planar wrt {7t,tt'); 
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if a is planar wrt (vr, vr') 



a 



m=l 
X 



m=l 



(5-8) 

with the decompositions of n^^a and n'^^a into their products of cycles : 

vr^V = 7ri7r2 . . . vr^ , vr'" V = 7fi7f2 . . . vf „ (5-9) 

sitc/i that: 

Remark 5.2 From the definition, one can see that o"(rm,a„) = vr(rm, 1) and cr{rrn,am) = 
7r'{fm, 1) for any m. 

5.1 Properties of the C^'s. 

Lemma 5.1 The matrices are symmetric. 

proof: 

It comes directly from the definition. • 



Lemma 5.2 



proof: 

One has: 



T.Ct = Id (5-11) 



Wl^,{xuyu X2, 2/2, ... , Xk, yk) = E (^a.-x IT ^^i^i^ y-ii)) (5-12) 



Let us shift aU the x's by a translation a and send a —>■ oo, i.e. replace all the 



Xj's by Xi + a. In the LHS, only the 5-terms of Definition |5.1| survive in the limit 
a — * oo, and thus the LHS tends towards the identity matrix. In the RHS, notice 
that Wi{xi + a,y„(^i)) 1. And C^, which depends only on the differences between 
independent of a. • 
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5.2 Some commutation properties 

Definition 5.3 Let A4''{x,y,C,,f]) be the k\ x k\ matrix defined by: 

Let A^''\xi,yi, . . . ,Xk,yk) be the k\ x k\ matrix defined by: 

{Al^l{xi,yi, ...,Xk, yk) ■= J2i Xiy^{i) 
A'-^yixi, yi,..., Xk, yk) ■= 1 if ti^tt'^^ = transposition 
A'-^yixi, yi,..., Xk, yk) ■= otherwise 

Theorem 5.2 



(5-13) 



(5-14) 







[M\x,y,^,v),Ct{x,y)]=0 


and 










[M'{x,y,^,r]),W\x,y)]=0 


proof: 

Let us define: 


M{x, 


y,^,v) ■■= M{Nx,y,N^,r]) 



(5-15) 



(5-16) 



(5-17) 



and 



Pm 



Wty{xuyi,...,Xk,yk) := ( H \^Pm,i+ Tr H 



\m=l 



\N{M,-x,^JiM,-y,^^^ 



It was proven in [15] that: 

[M'ix,y,^,rj),W'ix,y)]=0 



(5-18) 
(5-19) 



Now, in the large N limit, the factorization property [5] < Tr Tr >~< Tr >< 
Tr >, implies: 

yV^y{xi,yi,X2,y2, ■ ■ ■,Xk,yk) 

n 



JVnc.c,..(-'- -)-^y^l [x^^ y^^ 1/2, ... , Xfc, yk) 



(5 - 20) 



and using theorem 5J-, we have: 

>V^^,(xi,?/i,X2,|/2, . . ■,Xk,yk) 
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;V-c.c,ceK-i-)-fc 2/1, a;2, 2/2, • • • , x^, yt) fl (5-21) 



i=l 



Notice that 



(5 - 22) 



Thus: 



(5-23) 



Zl^^,7r,7r'(^a;i,2/i, A^a;2,2/2, . . . , Nxk,yk)Y[Wi{xi,y^(^i)) 

cr i=l 

Then, from [15], we have: 

k 

= [m\Nx, y, NC, v),CtiNxuyi, Nxk, yk)] fl W,{x„ y^^^^) (5-24) 



4 = 1 



In particular, choose a permutation a, and take the hmit where yi F(Xo--i(j)), you 
get in that hmit: 







M'iNx, Y{xl-i), v),C'ANxi, F(x.-i(i)), . . . , Nxk, Y{x^-i^k)))] (5-25) 



Since this equation holds for any potentials Vi and V2, it holds for any function 
Y{x), and thus the F(xj)'s can be chosen independentely of the Xj's, and thus, for 
any yi,...,yk,we have: 







M'^iNx, y, N^, r]),CaiNxi,yi, Nxk, yk) 



(5-26) 



Since it holds for any Xj's and ^, it also holds for Xi/N and C,/N. 



Corollary 5.1 



Va, [A^'\x,y),C'^ 



(5-27) 



proof: 

The corollary is obtained by taking the large ^ and rj limit of theorem ^ 
Appendix of [15]). • 



see 



30 



5.3 Examples: k = 2. 



(2)_/ ^1^11 1^22 ix,-x2){y^-y2) ] (^_2S) 



where Wij = Wi{xi,yj). 



C]^= i ^ {xi^x2){yi~y2) (5_29) 

V {xi-X2){yi-y2) ' 



f 2 _ {x^-x2){y2-y{) —I —Cl, 

^(12) — I 1 \—^^Id 



(5-30) 



6. Application: Gaussian case 



There is an example of special interest, in particular for its applications to string 
theory in the BMN limit [22], it is the Gaussian-complex matrix model case, V\ — 
V2 — 0. In that case one has E{x, y) — xy — 1, and thus: 

Wr{x,y) = ^^ (6-1) 
xy-1 

The loop equation defining recursively the W^s can be written: 
{xiyk - l)Wk{xi, yi,..., Xk, yk) = 

xWk-j+i{xj,yj, ...,Xk,yk) 

(6-2) 
Its solution is then: 

k 

Wk{xi,yi,. . . ,Xk,yk) = Cl''\xi,yi,. . . ,Xk,yk) ]J ^iML (6-3) 

CTGSfc i=l ^iy<T(i) ^ 

From the loop equation, one can see that Wk{xi,yi, . . . ,Xk,yk) may have poles 
only when Xi yj^ for any i and j. Because the C^-'s are rational functions of all 
their variables and because Wk has no singularity when Xi = Xj or y^ = yj, one can 
write: 

... / N NkiXi,yi,X2,y2:---:Xk,yk) 

Wk{xi,yi,...,Xk,yk) = pj-^ 7T (6-4) 

Y{i,j{xiyj - 1) 

where Nk{xi, yi,X2,y2,---,Xk,yk) is a polynomial in all its variables. 
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Moreover, the loop equation taken for the values Xk = or = shows that 
Wk{xi, yi,...,0, yk) = Wk{xi,yi, ...,Xk,0) = 0. Using the cychcity property of 
Wk{xi,yi, . . . ,Xk,yk), one can claim that it vanishes whenever one of its arguments 
is equal to 0. One can thus factorize the polynomial Nk{xi,yi,X2,y2, ■ ■ ■ ,Xk,yk) as 
follows: 

^ Qkixi,yi,X2,y2,---,Xk,yk)UiXiyi , . 

Wk{xi,yi,...,Xk,yk) = 7^-^ —r. (6-5) 

where Qki^i, yi, . . . , Xk, yk) is a polynomial of degree k — 2 with integer coefficient in 
all its variables. 

Notice that Qk{xi,yi, . . . , y~^i),yi, ... ,Xk,yk) = if a is not planar. 
As an example, we have: 

• for k = 2: 

W2{xi, yi,X2, 1/2) = „ ^^■^'^y^y'^ ^ Q2{.xi,yi,X2,y2) = 1 (6-6) 

Yli,Miyj - 1) 

• for /c = 3: 

W2.[xi,yi, X2, 1/2, 2:3, 1/3) = (2 - 2^ Xiyi+i + xiX2X2.yiy2y?,) -=r-( tT (6- 7) 

i Y{i,j[xiyj - 1) 

and 



Qz{.xi,yi, X2, 1/2, 3:3, yz) = (2 - Xiy2 - X2y-i - x^yi + xiX2X3yiy2y3) (6-8) 



7. Conclusion 

In this article, we have computed the generating functions of discs with all possible 
boundary conditions, i.e. the large N limit of all correlation functions of the formal 
2-matrix model. We have found that the 2k point correlation function can be written 
like the Bethe Ansatz for the 5-interacting bosons, i.e. a sum over permutations of 
product of 2-point functions. That formula is universal, it is independent of the 
potentials. 

An even more powerful approach consists in gathering all possible 2k point corre- 
lation functions in a /c! x k\ matrix W^. We have found that this matrix satisfies 
commutation relations with a family of matrices Ai'' which depend on two spectral 
parameters, and are related to the representations of U{n) [15]. We claim that the 
theorem p.2| is almost equivalent to the loop equations, and allows to determine W''. 

It remains to understand how all these matrices and coefficients C^- are related 
to usual formulations of integrability, i.e. how to write these in terms of Yang Baxter 
equations. For instance, the similarity with equations found in Razumov-Stroganov 
conjecture's proof [6] is to be understood. 
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One could also hope to find a direct proof of theorem |3]^, without having to 
solve the loop equations. In other words, we have found that the 2/c-point function 
can be written only in terms of Wi, while, in the derivation, we use the one point 
functions Y{x) and X{y) although they don't appear in the final result. 

The next step, is to be able to compute the expansion of those correlation 
functions, as well as the large limit of connected correlation functions. We are 
already working that out, by mixing the approach presented in the present article 
and the Feynman graph approach of [9] generalized to the 2-matrix model in [10]. 

Another prospect is to go to the critical limit, i.e. where we describe generating 
functions for continuous surfaces with conformal invariance, and interpret this as 
boundary conformal field theory [21]. 
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Appendix A Practical computation of f^j 

In this appendix, we build a set of trees in bijection with and use it in order to 



compute practicaly the weights fa defined in Definition 



Definition A.l Let Tk he the set of trees defined as follows: A tree T belongs to 
T E Tk, and is called a k-planar tree if and only if: 

• its root is labelled k + 1; 

• it is composed of k+1 vertices, labeled by [1, . . . , A; + 1]; 

• it has k edges which can be either upgoing or downgoing; 

• its vertices have valence 1,2 or 3, and are of one of the following eight possi- 
bilities, in which the point m denotes the origin of the branch containing i: 

— two trivalent vertices : 

^j+l<mm 

(1-1) 

— four bivalent vertices : 

^i+l<m m 

, "^i^ (1-2) 



:i-3) 



m i-l>m 

m m-1 ^ i>m 

two monovalent vertices corresponding to the leaves of the tree: 

i<m m 



i>m 



:i-4) 



Remark A.l Those trees are often called planar binary skeleton trees. 



Remark A. 2 One can see that the first edge is necessarily upgoing, and its extremity is 
necessarily 1. 



k+1 



Theorem A.l There is a bijection between Tk o-nd the set of planar permutations 
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proof: 

We build explicitly this bijection between and Tk-i. 



Consider a planar permutation a G ^k-i- Planarity means that a defines a 
partition of the disc into faces of two kinds. Let us say that faces which correspond 
to cycles oi Sa are colored in white, faces which correspond to cycles of a are colored 
in black. 

Decompose a and Sa into products of irreducible cycles: 

a = aia2...ai , Sa = aia2. ■ .a^ (1-6) 

and we assume that ai and ai contain xi. 

Because of planarity, we can define a distance of faces (i.e. cycles) to the face 
ai, as the number of edges one has to cross for going from a face a^ or ai to ai, and 
call it -D(o"j) or D{ai). 

We also define the "origin" of a face, noted m{ai) or m{ai), as follows: If the 
face is ai, we define m{ai) = k, otherwise, because of planarity, there is only one 
neighbouring face which is at smaller distance of ai. Because of planarity, those two 
faces share at most one x, and the origin is defined as the label of that x. 

Thus, each face has a color, white or black, a distance D, and an origin m. 

Now, to every face we associate a branch as follows: 
• to a white face, aj, i.e. a cycle of Sa, noted 

= ^i,2, • • • , ijj.) , = m{aj) , a{ij,n) = ij,n+i - 1 (1-7) 

we associate the upgoing branch ij^i — > ij^2 ■ ■ ■ —>■ ijj^ 




(if Ij = 1, the sequence contains only one vertex ij^i = m{aj) and no edge). 
• to a black face, aj, i.e. a cycle of a, noted 

we associate the downgoing branch ij^i —>■ ij^2 —>■...—>■ 




(1-10) 
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m{aj) and no edge). 
cr{ii,n) = ii,n+i - 1 (1-11) 

(1-12) 

By definition of the origin m of a face at distance D, the origin of a branch is 
necessarily a vertex on a branch at distance D ^ 1, and from planarity, it cannot be 
a vertex on any other branch. Thus, there is a unique way to attach all branches to 
their origin, and we obtain a tree, which belongs to 

Inverse bijection: 

On the other hand, let us consider a A; — 1-tree. One can build a permutation 
a e as follows: the image of an element of (1, . . . , /c — 1) is : 

• its descendant along a downgoing propagator if it exists; 

• the origin of the downgoing branch to which it belongs in the other cases. 

Because of the form of the vertices, the upgoing branches are necessarily the 
cycles of Sa. And since the branches form a tree, it implies that two faces touch one 
another through zero or one edge. Thus the permutation a is planar. 

It is easy to see that this application is the inverse of the preceding one. • 

Example: 

Let us carry out explicitly step by step this building for the permutation a G S12 
introduced earher. Notice that it is enrooted in 12 + 1 = 13. 




(1-13) 



(if Ij = I, the sequence contains only one vertex ij^i = 
• to the first face ai, 

5-1 = (ii,i,ii,2,...,h,rj , ii,i = i 

we associate the upgoing branch k — > ii^i — > ii^2 — > . . . 
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Considering the last non trivial cycle ag = (10, 11), one obtains finally the tree 
corresponding to (|4-18|) : 



/ \X \^ (1-14) 

13 3, 



Corollary A.l = Cat {k) , where Cat [k) is the k'th Catalan number. 

From now, one can see how to simply compute the weight f„ associated to a 
permutation a. 

Consider a planar permutation a G and its representation under the form of 
a tree T G #7^. Associate a weight to every edge of the tree as follows: 

• To every downgoing edge of a cycle of a, one associates the weight 

"11,1+ 1 

9in,l,in,j + l,in,j + 2 ) 

• To every upgoing edge of a cycle (T„ of Sa, one associates the weight 

■ii.j 

Then, the fa- is the product of all the weights of edges composing T. 
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